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ABSTRACT 

A correlation between the surface density of foreground galaxies tracing the Large Scale Struc- 
' ture and the position of distant galaxies and quasars is expected from the magnification bias 

Q\ . effect (Canizares 1981). We show that this foreground-background correlation function wjb 

can be used as a straightforward and almost model-free method to measure the cosmological 
Q_i' ratio Q/b. For samples with appropriate redshift distributions, Wfb oc $7 < 5g >, where 5 and g 

are respectively the foreground dark matter and galaxy surface density fluctuations. Therefore, 
Q/b oc Wfb/w, where w =< gg > is the foreground galaxy angular two-point correlation func- 
tion, b is the biasing factor, and the proportionality factor is independent of the dark matter 
power spectrum. Simple estimations show that the application of this method to the galaxy 
and quasar samples generated by the upcoming Sloan Sky Digital Survey will achieve a highly 
(~\| ' accurate and well resolved measurement of the ratio Q/b. 

^ Subject headings: Cosmology; gravitational lensing; biasing factor 

o\ ■ 

O ! 1. INTRODUCTION 

G\ ■ 

Q\ . Weak lensing promises to be one of the most effective ways to study the properties of Large Scale Structure 
*^ | (LSS) in the next years (Kaiser & Squires 1993; Kaiser 1998; Van Waerbeke, Bernardeau & Mellier 1999; 
O ■ Bartelmann & Schneider 1999). Most of the weak lensing methods proposed so far rest heavily on the analysis of 
JL \ background galaxy distorsions (Kaiser 1992; Bernardeau et al. 1997, Schneider et al. 1997), as they potentially 
5_i ■ contain finer detail about the LSS than the fluctuations of the background number counts, which are affected by 
\ intrinsic clustering. However, detecting and analyzing the typical shear expected from the LSS offers technical 
j3 ■ problems difficult to underestimate (Kaiser, Squires & Broadhurst 1995; Bonnet & Mellier 1995; Van Waerbeke 
J> , et al. 1997; Kaiser 1999; Kuijken 1999; although see Schneider et al 1998). In addition, a number of surveys 
will be available in the near future, as the Sloan digital Sky survey ( punn fc Weinberg 1995| ), NOAO (Januzzi 
et al. 1999) etc. with vast amounts of data which may not be optimal for shear analysis due to the imaging 
| pixel size or typical seeing. It is therefore necessary to develop methods which are able to tap the wealth of 
cosmological information contained in such surveys without relying exclusively on image distorsion analysis. 
Such an approach is provided by the background-foreground correlation function Wfb ( |Bartclniann fc Schnei- 



> 



X 



dcr 1991 ; Bartelmann fc Schneider 1993| ). The value of this statistic has been calculated by several authors 



using linear and nonlinear evolution models for the power spectrum evolution (Bartelmann 1995; Villumsen 
1996; Sanz, Martinez-Gonzalez & Bemtez 1997; Dolag & Bartelmann 1997; Moessner & Jain 1998; Moessner, 
Jain & Villumsen 1998). The two most obvious cases in which it is possible to measure Wfb are galaxy-galaxy 
correlations and quasar-galaxy correlations. The detection of the latter has a long and controversial history 
(for a discussion see Schneider, Ehlers & Falco 1992, Bemtez 1997), but it seems to be already well established. 
However, due to the scarcity of complete, well defined quasar catalogs not affected by observational biases, the 
results have low signal-to-noise and are difficult to interpret (Bemtez, Sanz & Martinez-Gonzalez 1999). The 
value of the expected amplitude for the low-z galaxy-high-z galaxy cross-correlation is rather small, and hard 
to measure within typical single CCD fields (|Villumsen, Freudling fc Da Costa 1997 ). Only with the advent of 



deep, multicolor galaxy samples and reliable photometric redshift techniques it has been possible to detect this 
effect (Herranz et al. 1999). 

To interpret the measurements of Wfb using the calculations mentioned above, it is necessary to assume 
a certain shape for the power spectrum. Unfortunately it is still far from clear whether the most popular 
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ansatz, that of Peacock and Dodds (1996) — or any other for that matter — provides an accurate fit to the LSS 
distribution (see e.g. Jenkins et al. 1998). It is thus desirable to develop methods whose application is not 
hindered by this uncertainty. An example is the statistic R (Van Waerbeke, 1998), which combines shear and 
number counts information, and can be used to measure the scale dependence of the bias. The value of R 
is almost independent of the shape of the power spectrum if the foreground galaxy distribution has a narrow 
redshift range. Here we show that something similar can be achieved with Wfb, with the additional advantage 
of being able to do without the shear information in those cases where the latter is difficult to obtain. 

It follows from the magnification bias effect (Canizares 1981, Narayan 1989, Broadhurst, Taylor & Peacock 
1995) that the surface number density of background population rib is changed by the magnification \i, associated 
with a foreground galaxy population with number density n as gb = (a — l)Sfi, where gb is the perturbation in 
the galaxy surface density rib {,9b = n b/ < rib > —1), ol is the logarithmic slope of the number counts and in 
the weak lensing regime fi ~ 1 + 8ii , S/j, « 1. Since fi ~ 1 + 2k and k, the convergence, is proportional to the 
projected matter surface density S, it follows that gb oc 5k oc 5T, oc 0,5, where 5 is the dark matter surface density 
perturbation. Therefore Wfb =< ggb >oc Q < g5 > and assuming linear and deterministic bias Wfb oc £lb~ l w 
where w is the two-point galaxy correlation function for the background populations, and the biasing factor b is 
defined by w = b 2 wss- This result provides a straightforward, virtually model independent method to estimate 
the ratio Vt/b. Williams & Irwin 1998 arrived to a similar expression using a phenomenological approach. The 
redshift distorsion method provides the quantity (3 oc fl os /b (Dekel 1994), so combining both one can estimate 
f2 and b separately. 

The outline of the paper is the following. In Sec. 2 we show rigorously that under reasonable assumptions 
Wfb oc Qb~ l w. Sec. 3 explores the application of this method to future and ongoing surveys and Sec. 4 
summarizes our main results and conclusions. 

2. FOREGROUND-BACKGROUND CORRELATIONS 

Let us consider two populations of sources: a background one (e.g. quasars or galaxies) and a foreground one 
(e.g. galaxies), placed at different distances A with p.d.f.'s Rb{X) and R(X), respectively. In Sanz et al. (1997), 
the background- foreground correlation Wfb was calculated as a functional of the power spectrum P(X,k) of 
the matter fluctuations at any time (A is the comoving distance from the observer to an object at redshift z, 

A = (i+ffl^+n andn<l,A = 0): 

w fb (6) = ( ab -l)C^(0), (1) 



C^s(6) = 120^ d \b(\,s0)T b (\)R(\)—^—T(\,se), (2) 

T(\, s e)^±-J™dkkP(\,k)j (k S 9), 1 _ {l X _ n)x2 , (3) 

where 6(A, sO) is the bias factor for the foreground galaxies (assumed to be linear, non-stochastic but possibly 
redshift and scale dependent), is the slope of the background source number counts and C^s{0) = 2 < 

/j,(4>)5((f> + 9) > is the correlation between the magnification and the mass density fluctuation. T&(A) is the 
lensing window function given by equation (7) in Sanz et al. (1997): 

The angular two point correlation function for the foreground population w can be obtained as 

w (9)= f 1 d\b 2 {X,s9)R 2 (X)[l - (1 - n)X 2 ]T(X,s9), (4) 
Jo 

where w(9) =< g((p)g((f> + 9) > is the foreground angular galaxy-galaxy correlation function. 

If we assume that the foreground and background galaxies are concentrated at the 'effective' distances Xf and 
Xb (a good approximation for realistic, nonoverlapping redshift distributions, Sanz et al. 1997), the previous 
formulas can be rewritten as 

w fb (X f , X b , 9) ~ K - l)12fi^ r r(A / , s f 9) (5) 
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w(X f , AX f , 9) ~ b 2 ^-[l - (1 - 0)A/]r(A / , (6) 
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where X C (A&, A/) = ofD Jh * s ^h e critical surface mass distance defined by the angular distance D, aj = izn^j^s 

is the scale factor, and AAy ~ l/Rf(Xf) is the width of the p.d.f. Rf(X) (e.g. the FWHM for a Gaussian 
distribution). Dividing Eq. (5) by Eq. (6) and simplifying one obtains 



w b 
i.e. the two correlations are proportional. The proportionality factor Q has the form 



-QAz f (a b -l)- (7) 
w b 



n = QX f (i-x f /x b )[i-(i-n)x f x b ](i-x f ) 

[l-(l-0)A / ][l-(l-fi)A^ W 

Fig. 1 (top) shows the dependence of Q on z b and Zf for the = 0.3, A = case. For the J) + A = 1,A>0 
case, there is not an explicit form for the angular distances, and the corresponding expression for Q is (Fig. 1, 
bottom) 

Q _ 12(1 + , f) 



X c (z f ,z b )Ji + Qz f + A[(l + Zf y 2 - 1] 



Note in Fig. 1 that for Zf < 0.2 and z b > 1 the value of Q in both open and flat geometries changes very 
slowly with z b and almost linearly with Zf. To better understand this behavior from Eq. @ let's assume that 
Zf << 1 < z b . In that case, A/ ~ sj ~ gZ/, a/ ~ 1 — z/ and Q ~ Therefore, Wf b is almost independent of 
X b (Fig. 1) and roughly Wf b ~ 3zjAzj-(ab — l)r26 _1 u; . Fig. 2 also shows that for the same redshift range, the 
amplitude of Q is also approximately independent of £1, allowing an empirical, model-independent estimation 
of fl/b. 

3. PRACTICAL APPLICATION 

The main quantity which determines the signal in the measurement of the angular correlation function is the 
excess (or defect) in the expected number of galaxy pairs. For a bin with surface A b i n at distance 6 this number 
will be 

AN {6) « (±)N f n b A bin (9)w fb (9) (10) 
Therefore, the signal-to-noise of the detection will be roughly 

S AN 



N JN 



N f ^n b A bin w fb {9) (11) 



We have not included the scatter due to the clustering of foreground galaxies and background galaxies, but it 
is unlikely that this effect will increase the noise over Poisson more than a factor of a few for any realistic case. 
The quantity n b A b i n is the number of background galaxies per bin. If we assume radial concentric bins of width 
A# and Wf b = Aj b 9~^, the signal to noise in each bin will be 

S 



-~A} b ^2irN f n b A06°- b -i (12) 

Since 7 is typically 0.7 — 0.8, the efficiency of the method will decrease very slowly with radius, allowing to map 

0/6 up to very large scales. 

The SDSS will obtain ~ 10 6 spectra for a population of < z >~ 0.1 galaxies ( |Gunn fc Weinberg 1995 ) 



forming a splendid foreground sample. It may be assumed that its projected angular correlation function will 
be similar to that of the APM catalog, which has an amplitude A = 0.44 at 1' and a slope 7 = 0.668 flMaddox 



et al. 1990D - The SDSS will also obtain S/N ~ 10 u'g'r'i'z' photometry for 10 8 galaxies in a 10 4 degree region, 



for which photometric redshifts will be estimated. The background sample can be formed by those galaxies 
with z > 0.4 and < z >~ 0.7 — 0.8, with n b ~ 1 — 1.5/d 7 . The value of Q cannot be considered constant for 
these two samples (see Fig. 1) because of the range of redshifts involved. Luckily, the existence of spectroscopic 
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redshifts for the foreground sample will allow to select subsamples with an extremely thin redshift distribution 
for which Q is approximately constant, each yielding independent estimates of Q/b which can afterwards be 
combined together. A rough estimate of the expected result using Eq. (12) gives 



SDSS 




(13) 



(we have assumed Nf = 10 6 , n& = 1.5, Azf = 0.05, Q = 0.25, aj> — 1 = 0.5) This shows that it will be possible to 
map 0,/b with an extremely good combination of resolution and accuracy using the SDSS. At large radii, the 
bias factor is expected to become constant ( Poles 1993|) , which means that the bin size can be made as large 
as desired without being affected by the scale dependence of b and therefore substantially decrease the error in 
the estimation of 0,/b. It is obvious from the above numbers that the main source of errors will not be the shot 
noise, but contamination due to low redshift galaxies which may sneak into the high-redshift sample creating 
a spurious correlation. This contamination can be very effectively minimized by applying a Bayesian threshold 
(Bemtez 1999). In addition, it will be possible to accurately quantify and correct this effect using a calibration 
sample with spectroscopic redshifts. 

The SDSS also plans to obtain spectra for a sample of 10 5 red luminous galaxies with < z >~ 0.4 and 10 5 
QSOs. Due to the much lower density of the background sample, it will not be possible to attain such a good 
spatial resolution in the measurement of Q/b, but these samples will allow to trace the evolution of biasing with 
redshift up to z ~ 1. However, experience shows that in this case the utmost attention will have to be paid 
to eliminating or accounting for the observational biases in the QSO detection and identification procedure, 
which can totally distort the estimation of w^f (Ferreras, Bemtez & Martinez-Gonzalez 1997, Bemtez, Sanz & 
Martinez- Gonzalez 1999). 



4. CONCLUSIONS 

A correlation (positive or negative) between the surface density of foreground galaxies tracing the Large 
Scale Structure and the position of background galaxies and quasars is expected from the magnification bias 
effect (Canizares 1981). We show that this foreground-background correlation function wjb can be used as a 
straightforward and almost model-free method to measure the cosmological ratio £l/b. 

For samples with appropriate redshift distributions, Wfb oc Q < Sg >, where 5 and g are respectively the 
foreground dark matter and galaxy surface density fluctuations. Therefore, Q/b oc Wfb/w, where w =< gg > 
is the foreground galaxy angular two-point correlation function, b is the biasing factor, and the proportionality 
factor is independent of the dark matter power spectrum. 

Simple estimations show that the application of this method to the galaxy and quasar samples generated by 
the upcoming Sloan Sky Digital Survey will achieve a highly accurate and well resolved measurement of the 
ratio Q/b. 
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Fig. 2. — Dependence of the factor Q (Eqs. 
and Zb. The thick lines correspond to Q, + A : 



;. g, 11) on Q for different redshifts of the foreground and background populations Zf 
— 1 universes and the thin ones to A = 0. 



